function a = gen_a(M, Delta, theta)
%GEN_A Summary of this function goes here
%   Detailed explanation goes here

vect = 0:M-1;

a = exp(i*2*pi.*vect'*Delta*sin(theta));

end

